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Abstract 

OO 

The entanglement characteristics of two qubits are encoded in the 
invariants of the adjoint action of SU(2)(g)SU(2) group on the space of 
density matrices ^3_|_ , defined as the space of 4 x 4 positive semi-definite 



^) Hermitian matrices. The corresponding ring C[*P+] su ( 2 )® SU(2) Q f 

polynomial invariants is studied. The special integrity basis for the 
ring C[^p+] su ( 2 )® SU(2) ^ - g described anc [ constraints on its elements 
due to the positive semi-definiteness of density matrices are given ex- 

S^ plicitly in the form of polynomial inequalities. The suggested basis 

5-j is characterized by the property that only a minimal number of in- 

variants, namely two primary invariants of degree 2, 3 and one sec- 
ondary invariant of degree 4 appearing in the Hironaka decomposition 
of C[*}3+ jSU(2)®SU(2) _ are g U bj e ct to the polynomial inequalities. 



1 Introduction 

According to the quantum theory non-locality of a quantum word manifests 
itself in a way that is very different from the intuitive classical views. At 
the very outset of quantum epoch reflections on that fact create a variety 
of paradoxes beginning from the Eistein-Podolski- Rosen paradox and the fa- 
mous neither dead nor alive Schrodinger cat [U El E]- Only towards the end 
of the XX-th century after the advances of technology, when a manipula- 
tion with quantum coherency became reality, the pragmatic approach to the 
problem rises the question of a practical usage of quantum non-locality Time 
for the realization of quantum communications and creation of the quantum 
computer has to come [I]. 

The difference between quantum and classical correlations has a very 
transparent mathematical background. One can already see its roots com- 
paring the basic states of classical and quantum computers; bits and qubits. 
While an arbitrary n-bit string can be transformed into other by the so-called 
"local transformation" , acting on its constituent bits, in quantum case this is 
true for one qubit states only. In other words the action of "local transforma- 
tion" cease to be transitive for multiqubit systems [51 E]- The action of local 
transformations splits the space of an arbitrary quantum system into the 
equivalence classes in way that each class is characterized by a different non- 
local properties [7j. Therefore the problem of classification of non-localities 
in a system of n-qubits reduces to the mathematical problem of description 
of orbits of "local" group action on the space of states [H [9]. The corre- 
sponding orbit space, £ n , is termed as "entanglement space" [51 E]- For its 
characterization the mathematical formalism based on the classical theory of 
invariants (c.f. [101 E]) is often applied. In this approach, in order to sepa- 
rate orbits, i.e., to introduce coordinates on S n , the polynomials in elements 
of the density matrices, which are invariant under the local transformation, 
are used. 

The entanglement space has highly nontrivial geometric and topological 
structure [HI [12] . Especially, a complexity of £ n is steeply rising with number 
of qubits growing up. This makes computations very tedious. However, for 
the lowest, 2-qubits system, the approach based on the classical theory of 
invariants allows to obtain a series of important algebraic results, clarifying 
the properties of S 2 EH E]. 

There is one further complication with the description of £ n . In virtue of 
a physical requirement the density matrices should be positive semi-definite 



[T5"l [TBI [T7| . Therefore, the space of local group action is not a linear space, 
but represents a certain semi-algebraic variety *P + . This circumstance should 
be taken into account applying the classical theory of invariants for con- 
struction of the orbit space. In the present article this problem is analyzed 
and detailed solution is given for the case of 2-qubits. With this aim the 
semi-defmiteness of density matrices is formulated explicitly in the form 
of polynomial inequalities in the scalars of the adjoint action of the group 
SU(2) ® SU(2) . Apart from this, the integrity basis for the polynomial ring 
C[*}3 + ] su ( 2 ) (g,su< - 2 ) , that includes the minimal number of elements subject to 
the above inequalities, will be presented. 

Our plan is as follows. We start, in paragraphs 2 and 3, with the brief 
review of necessary notions from quantum mechanics and put them into 
context suitable for the characterization of entanglement within the classical 
theory of invariants. Further, in the 4-th paragraph the system of polynomial 
inequalities in the Casimir operators of enveloping algebra of SU(4) , that 
describes the space *}3 + is derived. Regarding to these inequalities, in the last 
paragraph the integrity basis for the ring C[^3 + ] SU( - 2 ^ su ^ 2 - ) is constructed. 

2 The space of states 

A generic mixed state of n-level quantum system is described by n x n complex 
matrix, the density matrix g [151 [16], satisfying following conditionsj 

i. Hermicity - g = g + , 

ii. finite trace - Ti(g) = 1 , 

iii. positive semi-definiteness - g > . 

Mixed states form the subspace *P + , of the space of Hermitian nxn matrices. 
It is instructive, before considering a generic n-level system, to start with the 
simplest two-level quantum mechanical model. 



1 The special class of idempotent matrices, g 2 — p , corresponds to the so-called pure 
states, whose description reduces to the usage of rays in a Hilbert space. A mixed state is 
a mixture of pure states. 



2.1 Qubit 

In quantum theory of information an abstract quantum mechanical model 
with two levels takes a special place and independently of its physical real- 
ization carries the universal name - qubit. 

The qubit state is given by a density matrix that coincides with the 
standard density matrix of the non-relativistic spin 1/2: 

e =-(l + a-a), (1) 

where a - set of the Pauli matrices and ex is defined as mathematical 
expectation: 

ex = Tr (erg) , 

In the representation (fTl) requirements (i.) and (ii.) are taken into account 
by construction. The condition (iii.) restricts the parameter space of mixed 
states by a unit ball 

« 2 <1, (2) 

while for the pure states of qubit the expectation ex. lies on the Bloch 2-sphere 

a 2 = l. 



2.2 Qudit 

Analogously to the qubit the special terminology for d - level quantum sys- 
tem, a "qudit", has been introduced. The generalization of representation 
([I]) to the case of qudits reads [H] : 



'=s("^)Ft^«-*I. (3) 

where £ = (A) G R d _1 is d 2 — 1- dimensional Bloch vector. In the expansion 
([3]) components of the vector A = (Ai, A2, • • • , Ad2_i) represent the elements 
of the su(d) algebra normalized by conditions 

2 
AjAj = - Sijld + {dijk + i /yfejAjt, 



The explicit form of a-matrices is given below, in paragraph 5, formulaes (18 1 



Sij is the Kronecker symbol. The constants d^ and f^ are the so-called 
totally symmetric and antisymmetric structure constants of the algebra: 

1 i 

dabc = TTr({A a , Aft}A c ), f a bc = — TTr([A a , Af,]A c ), 

where 

{A a , Ab} = \ a \b + A(,A a , [A a , As] = A a A(, — AftA a . 

As for the case of a qubit, the properties (i.) and (ii.) of the qudit's den- 
sity matrix are already taken into account in the decomposition (J3J). The non- 
negativity requirement (iii.) imposes further, more subtle than pi), restric- 
tions. A complete characterization of qudit's Bloch vector space, B(R d _1 ), 
in an arbitrary dimension, is an open problem. However, some general prop- 
erties of this space is already known. Particularly, it can be shown that 
B(IR d _1 ) is a convex subset of a d 2 — 1-dimensional unit ball 

£ 2 < 1. 

It being know that all pure states are concentrated on its surface. More 
precisely, qudit's pure states are determined by the equation 



* a = i, £V£ = £ 



where 



(£ V £) k := \l ^LJl J_ dijk ^ 



d-2 



*3 



2.3 Composite states 



From the standpoint of quantum information theory it is mostly interesting 
to consider states composed of several qubits. According to the quantum 
theory axiom on composite systems jl] , the space of states of system, which 
is obtained by joining two systems A and B , represents a subspace of the 
tensor product of their individual Hilbert spaces Ha and Hb '■ 

ncH A ^n B . (4) 

The definition (El) in conjunction with the superposition principle is a 
source of an existence of correlations in the joint system, which do not have 
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5> ifi £<g>ef, 


w j > , 


Z)^ 


3=1 




3=1 



any classical analog. If a mixed state g, describing the joint A + B system, 
admit (not necessary in unique way) representation of the form 

1, (5) 

where gf and gf are density matrices of subsystems, then this joint state 
is called separable [7]. For such a state correlations between subsystems are 
classically conceivable. But the states (p| are far to exhaust all possible states 
of combined system. The states that can not be written as (pi), are called 
entangled. 

For a pair of r and s-qudits it is useful to represent the density matrix in 
the so-called Fano form 119, 201: 



/ r z -l s^-l H-ls'-l \ 

g = — I I ra + ^2 a i ^ ® ^ + Z bi Ir ® Ti + Z Z Q -^ ® r -? ) • ( 6 ) 

\ j=l i=l i=l j=l / 

In ([6]) matrices A, and Tj are basis elements of the su(r) and su(s) algebras 

respectively. The real (r 2 — 1) x (s 2 — 1) matrix C = \\cij\\ is so-called 

"correlation matrix". Meaning of parameters a = (oi, . .. ,a r 2_i) and b = 

(pi, . . . , 6 S 2_ 1 ) becomes clear after performing the partial trace operation [5]: 

gM := Tr B (g) = -(I r + a ■ A) , ^ := Tr^) = -(l s +b-r). 

r s 

The vectors aab are Bloch vectors for subsystems whose states are describing 
by matrices g^ and g^ B \ respectively. 

The entanglement properties of density matrices d6J) as well as more gen- 
eral multipartite systems admit formulation in terms of invariants of the 
so-called local groups [9]. In the next paragraph the corresponding notions 
will be introduced. 



3 The entanglement space 

3.1 The local invar iance 

On the space of density matrices of n-level system the group SU(n) acts in 
adjoint manner 

Q -> Q ' = U ] gU. (7) 

6 



If a quantum system is obtained by combining r-subsystems with m, n 2 , . . . , n r 
levels, the non-local properties of the resulting composite system can be put 
into a correspondence with a certain decomposition of the unitary operations 
in (|7j). Namely, from all unitary actions we separate the group of so-called 
local unitary transformations (LUT) 

SU( ni )®SU(n 2 )®---®SU(n r ), (8) 

acting independently on the density matrices of each subsystems 

q^ -)• Q {ril) ' = g ] Q {ni) g ^GSU(ni), z = l,2,...,r. 



Two states of composite system connected by the LUT transformations (18j) 
have the same non-local properties. The latter can be changed only by the 
rest of the unitary actions 



SU(n) 

n 



SU(ni) ® SU(n 2 ) ® • • • ® SU(n r ) 



n 

i=l 



generating the class of non-local transformations. 

As it was mentioned in the Introduction the action of LUT is not tran- 
sitive. The equivalence of states regarding the action (|8j) gives rise to a 
decomposition of the space of matrices into the equivalence classes (orbits). 
The union of these classes, i.e., the orbit space, is customary to call as the 
"entanglement space" £ n . 

3.2 Orbit space and local polynomial invariants 

The main motivation for studying of £ n is necessity to work out qualita- 
tive criteria and quantitative measures for characterization of non-locality in 
composite systems [SJ [12] . 

As it was mentioned above, a canonical method for description of the orbit 
space £ n is the theory of classical invariants [UJ . Within this approach, start- 
ing from the works by Linden and Popescu [9] , series of interesting results, 
which clarify the mathematical contents of the entanglement phenomenon, 
has been obtained. A considerable progress was achieved for pure states. As 
an example, we refer here to the construction of Hilbert series for a multipar- 
tite systems of qubits [21] and classification of pure entangled states based 
on the theory of hyperdeterminants [22] . 



Analysis of the orbit space for systems in mixed states is much more vague. 
The generic questions of construction of basis for rings of local invariants for 
mixed states have been considered in [131 [H]. With this aim the algorithmic 
methods of computer algebra were used (23] , [2A\ |j 

According to the theory of invariants [Uj , the ring of polynomial invari- 
ants C[V] G , of linear space V over the complex numbers C, under the action 
of a group G , represents the graded algebra 

fe=i 

where A^ is the space of homogeneous G-invariant polynomials of degree k. 

The special unitary groups SU(n) belong to the reductive algebraic groups. 
Their ring is finitely generated |llj . and C[V] G is Cohen-Macaulay type [26J. 
However, straightforward application of this construction to the problems of 
quantum entanglement is complicated by the fact that the space *$+ , on 
which the local group ^ acts is not a linear space. As it was already em- 
phasized in the Introduction, density matrices are the positive semi-definite 
and therefore the space of representations *P + is nonlinear semi-definite alge- 
braic manifold. Bellow we suggest the trick how to overcome this difficulty, 
exemplifying the problem in details for a system of pair of qubits. 

Let start with the construction of the ring C[%4 X 4] SU ^® SU( ' 2 ' ) of the ad- 
joint action invariants on the space of 4 x 4 Hermitian matrices 7^4x4- In order 
to define the ring C[ < P + ] su ^ 2 ^® su ^ 2 - ) , note that the space of positive-definite 
matrices *$+ is subspace of "H4 X 4 , which is invariant under the action of 
SU(4) . As we demonstrate below the subset *$+ admits representation via 
the set of polynomial inequalities [j 

Pa(£ 2 ,£ 3 ,£4)>0, a = 1,2, 3 (9) 

in three invariants, (£2, £3 and £4 , of the enveloping algebra of SU(4) group. 
From the other side, since £2, £3, £4 are at the same time invariants of 
SU(2) ® SU(2) , then it is possible to construct in C[^ 4x4 ] su(2)0SU(2) such 

3 Unfortunately, applications of the existing algorithmic methods, including the Grobner 
bases technique, to the analysis of the ring of polynomial invariants for multipartite systems 
is not effective due to the sharp growth of the number of algebraic operations with the 
increasing number of qubits. 

4 One can find a description of *$ + , similar to the given here, in [271 1281129] , 



a basis that includes these invariants. As result, having this basis and tak- 
ing into account the inequalities (|9]), we will be able to characterize the ring 
C[*}3 + ] su ( 2 )® su ( 2 - ) completely. According to the consideration given in the 
subsequent paragraphs a basis of the ring can be chosen in a way that only 
the primary invariants of degree 2, 3 and one secondary invariant of degree 
4 presented in the ring's Hironaka decomposition [TT] are constrained by the 
above polynomial inequalities ^. 

4 Non-negativity of density matrix 

To succeed in our program of construction of an optimal homogeneous basis 
for the ring C[*P + ] su ( 2 )® su ( 2 ) let us start with the discussion of positive semi- 
defmiteness of density matrices. Below the requirement of non-negativity will 
be formulated in the form of inequalities in invariants of the adjoint action 
of SU(n) group on *}3 + . 

4.1 ^)3 + in terms of Casimirs of SU(n) 

A Hermitian operator is positive semi-definite if and only if all its character- 
istic numbers are non-negative. The condition of non-negativity of Hermitian 
operator can be formulated solely in terms of coefficients of its characteristic 
equation: 

\I n x-g\ = x n - S lX n - 1 + S 2 x n - 2 - . . . + (-l) n S n = . (10) 



The coefficients Sk in ( 10 ) are given as the sums of principal minors of k-th 
order: 



v ii ... i k 

l<ii<...<«j.<n 



Sk= ^2 g 

l<ii<...<jfe<n 

Since the matrix g is Hermitian all its characteristic numbers x k are real 



n. 



roots of the characteristic equations (10). When x fc are positive then all Sk 
being the symmetric polynomials 

k 

Sk= Yi n^> 

l<£l<...<ifc<n j=l 

are non-negative real numbers. The inverse statement is true as well; the 
non-negativity of the coefficients S k provides the non-negativity of roots x k ■ 

9 



The Descartes theorem [30]: a number of positive roots (taking into account 
their multiplicity) equals to the number of signs changes in the sequence of 
the coefficients of the polynomial equation, gives proof of this observation 
(see e.g. [27]). 

So, non-negativity of density matrices can be written in the invariant way 
as condition of non-negativity of the coefficients of its characteristic equation: 

S k >0, k = l,...,n. (11) 

We give here, for further use, the explicit form of a few first coefficients 
Sk, written in terms of n-dimensional Bloch vector £ [27, 28J: 

1 n — 1 



1 (n-l)(n-2) 
3! n 7 



Ss = ii Kk T ; (l-3g-g + 2(j-Vfl-g), 



1 (n - l)(n - 2)(n - 3) . _ ^ ,* ^ * 

+ 3^(^^) 2 -6^|(^VO-(^VO). 
n — 6 n — 6 



Besides from the restrictions (jllj) there are upper bounds on Sk due to 
the normalization condition Tr(^) = 1 , Tr(£> fe ) < 1 , for k > 2 . Note, that 
the equality fulfils for pure states and the maximal values of Sf. are achieved 
for the equal eigenvalues Xi of density matrices. 

Finally, the positive semi-definiteness and normalizability conditions for 
density matrices of n-level system can be written as the following set of 
inequalities 

0< ? 1V I, 7 7— r, <1, fc = 2,...,n. (12) 

(n — lj(n — 2J . . . [n — k + 1) 

Coefficients Sk, k — 1, . . . , n of the characteristic equation are invariants 
under the adjoint action of SU(n) group. They are algebraically independent 
and can be represented via polynomials in the Casimir operators of the corre- 
sponding enveloping algebra. Below, the case n — 4 , related to the system of 



2-qubits, is considered in details and inequalities (12) are rewritten directly 



in terms of the Casimir operators of the enveloping su(4) algebra. 



10 



4.2 Restrictions on invariants of SU(4) 

The group SU(4) has three Casimir operators whose expressions in terms of 
the components of 15-dimensional Bloch vector £ (see the decomposition (J3J)) 
can be written as 

£ 2 = £-£, (13) 

£ 3 = £V£-£, (14) 

£ 4 = £V£-£V£. (15) 

Because for an arbitrary 4- level system the coefficients S 2 , S 3 and 64 of char- 
acteristic equation of density matrix are expressible via these Casimir oper- 
ators 

5 2 = l(l-£ 2 ), 

5 3 = ^-(l-3<Z 2 + 2<£ 3 ), 

16 

Sa = t^((1-3€ 2 ) 2 + 8€ 3 -12€ 4 ), 
zoo 



the set (12) reduces to the following constraints on SU(4) invariants 



< £ 2 < 1 , 
< 3£ 2 - 2£ 3 < 1 , (16) 

< (1 - 3£ 2 ) 2 + 8£ 3 - 12£ 4 < 1 • 



In the space spanned by invariants £ 2 , £3 and £4 inequalities (16) define 
the bounded domain depicted on the Figure [Tj 

5 The ring of local invariants C[^P + ] SU ( 2 )® SU ^ 

Consider the density matrix of two qubits parameterized in the Fano form: 

e = - [1 2 ®1 2 + a- a ®l 2 + l 2 ®b- a + djGi® aj] , (17) 

where 3-component vectors a = (ai,a 2 ,a 3 ) and b = (bi,b 2 ,b 3 ) are Bloch 
vectors of constituent qubits, o"j , i — 1,2,3 - the Pauli matrices making up 
the basis of algebra su(2): 

*1 = ( J I ) > ^ = ( ° ~ l J , CT3 = (J ^ J . (18) 

11 




Figure 1: The allowed region for the values of Casimirs €2, £3 and £4 
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The correlation matrix C of pair of qubits has 9 elements c^- i,j = 1, 2, 3 . 
At first, following our program of the construction of the ring of invariants 



C[*P jsu(2)<gisu(2) ou tii nec { j n paragraph 3.2, we identify the space of parame 



ters a , b and C with IR 15 , ignoring for a moment all restrictions due to the 
non-negativity of density matrices. Apart from this we linearize the adjoint 
action Q of the local group SU(2) <g> SU(2) 

V A -> V A = L AB V B A,B = 1,...,15, (19) 



with 15 x 15 matrix L G SU(2) ® SU(2) ® SU(2) <g> SU(2). 

So, our preliminary issue is to build up the ring of polynomial invariants 
of linear action of SU(2) ® SU(2) <S> SU(2) <S> SU(2) group on the linear space 
M 15 . Note that the linearization (19) allows to use the prompt from the 



Molien formula for the generating function of invariants for txq representa- 
tions of a compact group G [21] : 

M(q) = [ dfi G ——± 1| , (20) 

J G det||id-g7r G || 

where the integral is taken over the group G with the Haar measure &jjlq . 

The Molien function provides information on the polynomial invariants 
ring's structure. Firstly, its formal series in powers of parameter q, the so- 
called Hilbert-Poincare series: 



M(q) = Y J d k q k G Z[q] 



k>0 

points out the dimension, dk, of the space of homogeneous invariants of degree 



k. Secondarily, being a rational function, (20) admits (non-uniquely) for 
q < 1 , the representation 

spr degJ k 

M(q) - ^ k=0 q 



n^ = i(i-g degi ^) 



From this form of the Molien function on can conclude on the number and 
order of the primary Ki , i — 1, 2, . . . , n , and secondary Jj , % — 1, 2, . . . , r , 
invariants of the Cohen-Macaulay algebra 



C[Vf = Q)J k CiK u K 2 ,...,K v 



k=0 
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As computations show, the Molien function for mixed states of two qubits 
can be written as p31 [H] : 

1 + q A + g 5 + 3g 6 + 2q 7 + 2q 8 + 3g 9 + g 10 + q 11 + q 15 
{q)= (l-g)(l-g 2 ) 3 (l-g 3 ) 2 (l-g 4 ) 3 (l-g 6 ) ' ( j 



According to the result (21), a basis of ring consists from 10 primary invari- 



ants of degree 1,2,2,2,3,3,4,4,4,6 and 15 secondary invariants of degree 
4,5,6,6,6,7,7,8,8,9,9,9,10,11,15. 

More detailed information on invariants dependence on coefficients of the 



decomposition (17) can be extracted using the so-called method of many- 
parametric generating functions [21]. In our case the many-parametric gen- 
erating function depends not only on one parameter q, but is a function of 
three arguments, F(a,b,c). The contribution from variables a,b and Cy 
into the Molien function now is taken with the weights determined by each 
independent parameter, a , b and c respectively. 

It is worth to note that the generating function F(a, b, c) was found al- 
ready in the middle of 70-th of the last century [2U E2] , in connection with 
so-called problem of "missing index" which arose within issue of nuclei spec- 
trum classification. The corresponding mathematical formulation and solu- 
tion of the problem can be found, e.g., in [31]. Further, in our presentation, 
we will mainly follow the article [32j . 

Consider the space of all polynomials in fifteen variables Oj, foj and c^- 
i,j = 1,2,3. In virtue of the adjoint action of the local group the space of 
Bloch's parameters is decomposed into the irreducible representations of the 
SO(3) (g> SO(3) group. More precisely, the variables a*, 6j and Cy are trans- 
formed according to the representations T)\ x D , D x £) 1( and Di x D\ 
correspondingly. Since the subspace, P a> t,q[a>h h, %•], of homogeneous poly- 
nomials in variables ai,bi,Cij of degree s,t, q correspondingly, is invariant 
under the action SU(2) £§> SU(2) , all invariants C can be classified according 
to their degrees of homogeneity, i.e., C^ stq \ 

Consider, following the construction suggested in [32] , the set of invari- 
antsS 

3 invariants of second degree 

C^ = c ijCij , C^= ai a t , C^ = bA, (22) 



5 Below, everywhere it is assumed the summation over all repeated indices from one to 
three. 
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£(004) 


Cj Q Cj^C j a C j p , 


£i(202) 


— (li(ljCi { yCj ct , 


£(022) 


- babpCiaCifi , 


£(112) = 


tjj7j6 Qj g^,UjO Q ,Cj j gC^^, 



2 invariants of third degree 

4 invariants of fourth degree 

(24) 
(25) 
(26) 
(27) 

1 invariant of fifth degree 

C (113) = didaCjaCjpbp , (28) 

4 invariants of sixth degree 

C (123) = e ijk biC aj a a C0 k Cpibi , (29) 

C (204) = CLiCiaCjaCtfCkpak , (30) 

C ( ° 24) = biCriCajCftCfrkh , (31) 

C^ — £al3jO<aCl3ibiCjjCsjas , (32) 

2 invariants of seventh degree 

C = tijkbiCajdaCpkCpiC^iai , (33) 

C (124) = eapyaaCfobjOykCskCsibi , (34) 

2 invariants of eighth degree 

t- 7 = £ijkbiC a jC a lbiCfii t Cp m Cy m (l.y , (35) 

C = ^a/s-yaaCfSiCsiasc^kCgkCgibi , (36) 

2 invariants of ninth degree 

O ^aP^/^aC/IJiCSiClsC^jCgjCgkC^^Clij , \^'/ 

O = €ijk0iC a jC a i0iCfjkC/3 m C'ymC^/s( ) s j v^"J 

15 



From these invariants the basis of C[^J + ] su ( 2 )® su ( 2 ) can be build. As the 
criterion of its construction we choose the principle of usage of basis with 
the minimal number of elements involved in the definition of ^3+ . Having 
in mind this rule and noting that the space *}3 + is defined in terms of the 



Casimir operators of SU(4) group ( 13 )-(15 ), we expand £2,^3, £4 over the 
set of above introduced local invariants Q22|)-(24): 



e 2 = - (c( 200) + c(° 2 °) + c(° 02 )) , 

£4 = - \2(C {200) C {020) + C (202) + C (022) - C (112) ) + (C(° 02 )) 2 - £(° 04 )] . (41) 
6 



(39) 
(40) 



From equations (39)-(41) it follows that one can consider the Casimir op- 
erators £2, £3, £4 as the basis elements instead of scalars C'- 002 *' , C 1 - 003 ^ and 

£(112) _ 

Bear in mind this observation and using the results of [14J , where the ring 
C[R 16 jsu(2)®su(2) wag described, we define the following set consisting from 
10 primary invariants, including the Casimir operators €2, £3, 
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K 7 = C^ , 
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K 9 
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K w = C^ + C(° 24 ) , 






5 secondary invariants including the Casimir £4 
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Ji = U, 
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J 2 = c (113) , 






deg = 6 , 


j 3 = c( 204 > - c^ , 


J 8 = c^ , 
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J u = C^ , 




deg = 8 , 


J12 = c( 215 ) , 


Jl3 = C^ , 




deg = 9 , 


J4 = J1J2, 


Ju = C^ , 


Ju 


deg = 10 , 


J 5 = JlJ^i 






deg= 11, 


J 6 = ^2^3, 






deg = 15 , 


■h — J1J2J3 ■ 







£(020) 



(42) 



£(022) 



£(213) 



(43) 



£(036) 
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We conclude that the set of homogeneous invariants (42)-(43) represents 
the basis for the ring C[^] su ( 2 )® su ( 2 ) : 

cpp+] SU(2)BSU(2) = Q Jk C[Rh R2 ^ Kw] 
k=0 

under the condition, that two primary invariants K 2 , K 5 and one secondary 



invariant J\ satisfy the inequalities (16) 



6 Conclusion 

An essential issue of the quantum theory of information is qualitative and 
quantitative characterization of purely quantum correlations caused by the 
entanglement of quantum states. Theory of classical invariants provides tools 
for studies of the corresponding space of entanglement, i.e., the orbit space 
of action of the group of local transformations on the space of states of 
composite systems. For the case we are interesting in, system of two qubits 
in a mixed state, the local transformations of the density matrices form the 
SU(2) (g) SU(2) group. Its adjoint action, on the space of the Hermitian, unit 
trace matrices, identified with IR 15 , defines the principal orbit space 

o-= *! 

' SU(2)®SU(2)' 

with dimension 

dimO = 15-2 x 3 = 9. 

However, the orbit space defined in such a way is not the space of entangle- 
ment £2 ■ Due to the non-negativity of density matrices the space of physical 
states is *P + C M 15 . In the present article we suggest the description of *}3 + 
based on the polynomial inequalities in Casimir operators of the enveloping 
algebra su(4) . Furthermore, we show how these restrictions can be effectively 
taken into account constructing the basis for the ring C[*P + ] su ( 2 )® su ( 2 ) ? pro- 
vided for the Hironaka decomposition with only two primary invariants of 
degree 2, 3 and one secondary invariant of degree 4 constrained by the poly- 



nomial inequalities (16) 



Concluding it is important to emphasize that without the inequalities 
(16), the usage of local invariants for "coordinatization" of the space of en- 
tanglement £2 is not correct. We leave for a future publications analysis of 
those constraints consequences on the geometry of £2 C O . 
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